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A complete thermodynamical analysis of the 2 + 1 dimensional massless Gross-Neveu model is 
performed using the optimized perturbation theory. This is a non-perturbative method that allows 
us to go beyond the known large- N results already at lowest order. Our results, for a finite number 
of fermion species, N, show the existence of a tricritical point in the temperature and chemical 
potential phase diagram for discrete chiral phase transition allowing us to precisely to locate it. By 
studying the phase diagram in the pressure and inverse density plane, we also show the existence 
of a liquid-gas phase, which, so far, was unknown to exist in this model. Finally, we also derive N 
dependent analytical expressions for the order parameter, critical temperature and critical chemical 
potential. 

Q • PACS numbers: ll.10.Wx, 12.38.Cy, 11.15.Tk 
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I. INTRODUCTION 

The Gross-Neveu (GN) model [l[ has been extensively used as a prototype for quantum chromodynamics (QCD) 
and related issues. This is due to the fact that both models share some common features such as asymptotic freedom 
and chiral symmetry breaking (CSB). For this reason the GN model is useful as a toy model to test different tech- 
niques that can be ultimately used to tackle problems related to QCD phase transitions. At the same time, in the 
condensed matter physics domain, the two (1+1) dimensional GN model (GN2d) has been associated to polymers 
. including unidimensional molecules such as polyacetylene while the three (2+1) dimensional version (GN3d) 
' has been related to planar superconductors [J]. The applications concerning phase transitions within the GN model 
, are commonly carried out in the finite temperature and/or finite density domain where non-perturbative techniques 
must be employed. In these applications the most commonly used analytical non-perturbative technique is the 1/JV 
expansion Q where, for the GN model, N represents the number of fermionic species (see Ref. @ for a recent review). 
In general, this expansion is considered only at the leading order in what is known as the large- N approximation. 

In the large-TV approximation there are fundamental differences between the GN2d and GN3d models that are 
worth recalling. In dimensions d = 2 one observes chiral symmetry restoration (CSR) occurring via a phase transition 
of the second kind for high temperature (T) and small chemical potential (/i) values, while for low T and high fi the 
transition is of the first kind Q • One finds a tricritical point in the T — /i plane separating the second order transition 
line from the first order one, while metastable lines accompany the first order transition [3 Q . In the P—l/p plane (P 
being the pressure and p the density) one finds a phase diagram similar to the one generated by a Van der Waals liquid 
so that the CSB region corresponds to the "gas" phase while the CSR region corresponds to the "liquid" phase. It then 
follows that the first order transition allows for the appearance of a (mixed) liquid-gas phase. The chiral (dynamical) 
symmetry breaking happening in the massless GN2d with discrete symmetry at finite temperature, however, must 
be seen with care, since it is an artifact of the large- TV approximation. This is a consequence of well-known no-go 
theorems concerning that there should be no discrete symmetry breaking in one-space dimension Q . In this case the 
system's vacuum manifold allows for the appearance of kink-anti-kink configurations that are unsuppressed at any 
finite temperature (To| . The system becomes segmented into regions of alternating signs of the order parameter whose 
net average value becomes zero. At leading order, the l/N approximation misses this effect because the energy per 
kink goes to infinity asiV^ oo, while the contribution from the kinks has the form e~ N . The large- TV results for the 
GN model in 2+1 dimensions are rather different |4|]. First of all, as far a discrete chiral symmetry is concerned, the 
no-go theorem of one-space dimensions no longer applies (though now, in two-space dimensions, the no-go theorem 
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applies to the non-existence of a continuous broken symmetry at any finite temperature). The GN3d phase diagram 
produced by the large- TV approximation shows that the CSB/CSR transition is of the second kind everywhere except 
at T — 0, where it happens to be of the first kind. Within this approximation, there are no tricritical points lying in 
the T — n plane and no liquid-gas phase. The model behaves more like a planar superconductor with the transition 
CSB/CSR happening as a superconducting/normal one [4|]. Later, Kogut and Strouthos have used lattice Monte 
Carlo simulations to study the GN3d at finite TV They predicted that a tricritical point should exist at low 

finite values of T, but within the numerical precision of their simulations, they were unable to give its exact location 
(in earlier lattice simulations [T3| used to obtain the phase diagram for the GN3d at finite TV, some evidence for a 
tricritical point was also pointed out). No other attempts or approximations were able to improve on this situation. 
As far we are aware of, no evidence has been given so far for a possible chiral "liquid-gas" kind of phase. 

In an attempt to go beyond the simple large- TV approximation and in such a way that temperature and chemical 
potential effects could be considered in a simple approximation method, three of the present authors [l3| have recently 
considered the GN2d in the linear 6 expansion method (LDE), also known as the optimized perturbation (OPT). That 
study allowed for the inclusion of the first non trivial finite TV corrections to the complete phase diagram of the GN2d 
model. The main results of Ref. [l3| are the derivation of analytic non-perturbative expressions, containing finite 
TV corrections, scalar field expectation value , critical temperature T c (at p = 0), critical chemical potential /i c (at 
T = 0), as well as for the tricritical point (at T ^ and /t ^ 0). In the phase diagram, the predicted CSB region is 
reduced for finite values of TV. The OPT expression for T c predicts values that are lower than the ones predicted by 
the large- iV approximation which, in the light of the Coleman- Mermin- Wagner-Landau theorem [^] , can be viewed as 
an indication of convergence. 

Our recent success in treating the GN2d [l3| and the previous lattice Monte Carlo results on the GN3d, concerning 
the eventual existence of a tricritical point in the GN3d model [111, [l| , gave us the motivation to investigate the GN3d 
using the OPT method to fully study its thermodynamics in order to confirm, in an analytical way, the existence of a 
tricritical point. The OPT method is known for exactly reproducing the large- TV result for the effective potential (or 
free energy) already at the first non trivial order (ljjl [lj] . The perturbative computation of higher orders brings finite 
TV corrections and non-perturbative results are generated upon using a variational criterion. One advantage is that 
at any perturbative order one has complete control over the contributions, while the eventual technical difficulties are 
like the ones one should encounter in a traditional perturbative computation. The convergence properties of the OPT 
in critical problems associated to Bose- Einstein condensates have been proved [H, [l(| • It is worth mentioning that 
some of the most accurate numerical results regarding the critical temperature for weakly interacting homogeneous 
Bose gases have also been obtained with this method Concerning the GN3d model the results obtained in the 
present work include analytical equations for both T c (at /j, = 0) and fi c (at T = 0) as well as for the scalar field 
vacuum expectation value (vev) a c with finite TV corrections. Contrary to the GN2d case, these values appear to be 
higher than the predicted large-TV values. One of our most important results concerns the location of a tricritical 
point at finite values of T and p. Being able to specialize to any value of TV we choose TV = 3, which is the relevant 
value for QCD. In a preliminary work [181 ] . we already have shown that the unstable region in the T — p plane, which 
corresponds to the region inside the metastable lines that accompany the first order transition line, is rather small, 
thus explaining the difficulty in observing and locating the tricritical point in previous works. Here we extend that 
work and, from the (Landau) free energy, or effective potential, we obtain other relevant thermodynamical quantities 
such as the thermodynamical potential, pressure, density, etc. This allows us to obtain the phase diagram in the 
intuitively more accessible P — 1/p plane and that shows how important are the finite TV corrections in the GN3d. 
In fact, these corrections produce a phase diagram that is like a Van der Waals liquid and contrary to the large- TV 
predictions, we show that the model can display a mixed "liquid-gas" phase, which was previously unknown to exist. 

This work is organized as follows. In the next section we present the GN model. In Sec. Ill we present the OPT 
method and the interpolated GN model, evaluating the effective potential in this non-perturbative scheme. We show 
that, already at leading order, our results go beyond the known large- TV results. In Sec. IV we present the optimized 
results obtained from the effective potential at finite temperature and chemical potential. The dynamically generated 
fermion mass that we here associate with the auxiliar scalar field vacuum expectation value, the critical temperature 
and critical chemical potential for chiral symmetry restoration are evaluated and explicit analytical expressions for 
these quantities are obtained. In this same section we also discuss the complete phase diagram for the GN model in 
the T and p plane that then shows the presence of a tricritical point joining the lines of second order and first order 
chiral phase transitions, which we are able to locate precisely. In Sec. V we present other relevant thermodynamical 
quantities and show explicitly the existence of a mixed chiral symmetry restored/broken phase, the analogous of a 
liquid-gas phase in the P —1/p plane. The entropy, latent heat, and other important quantities are also evaluated. In 
Sec. VI we present the next order results, at T = and p — 0, that allows us to assess the convergence of the OPT in 
this model. Our conclusions are presented in Sec. VII. Three appendices are included to show some technical details 
and the renormalization for the interpolated model up to second order. 
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II. THE GROSS-NEVEU MODEL 

The Gross- Neveu model is described by the Lagrangian density for a fermion field tpk (k = 1, . . . , N) given by pj 

g 2 

where the summation over fermionic species is implicit in the above equation, with e.g. "ipk^k — Sfe=i Tpktpk- When 
nif = 0, the theory is invariant under the discrete transformation 1 

V> - TsV' , (2.2) 

displaying a discrete chiral symmetry (CS). For the studies of the model Eq. (|2.1j) in the large- N limit it is convenient 
to redefine the four- fermion interaction as g 2 N = A. Since g 2 vanishes like 1/N we study the theory in the large- iV 
limit with fixed A (see, e.g., |5|). 

At finite temperature and density the model can be studied in terms of the grand partition function given by 



Z(p, M ) = Tr exp [-(3 (H ~ M Q)] , (2.3) 

where (3 is the inverse of the temperature, /i is the chemical potential, H is the Hamiltonian corresponding to Eq. 
(12. ip and Q = J dx^klo^k is the conserved charge. Transforming Eq. (12. 3D to the form of a path integral in the 
imaginary-time (Euclidean) formalism of finite temperature field theory [19| . we then have 



N 

Z{fi,n)= / Y[D$ k Dij k exp{-S E $k,A}} , (2.4) 

^ k=l 



where the Euclidean action reads 



dr / dx 



ipk($ + M7o -m f )ip k 



2N 



(2.5) 



and the functional integration in Eq. (|2.4D is performed over the fermion fields satisfying the anti-periodic boundary 
condition in Euclidean time: ip k (x 1 r) ~ —ip k (x 1 T + (3). 



III. THE EFFECTIVE POTENTIAL FOR THE INTERPOLATED THEORY 



Let us now turn our attention to the implementation of the OPT method within the GN model. Usually, when 
employing this approximation one starts by performing a linear interpolation on the original model in terms of a 
ficticious parameter 5 (used only for bookkeeping purposes), which allows for further expansions. Accordin g to this 
OPT interpolation prescription [20j (for a long, but far from complete list of references on the method, see |2l|) the 
deformed four fermion theory reads [1 31 ] 

Csty, $) = $ k (i @) Vfe + (1 - 5) Ty^Vfc + ^(V^fc) 2 . (3.1) 

So, that at 5 — we have a theory of free fermions while at S = 1 the original theory is reproduced. Now, the 
introduction of an auxiliary scalar field a can be achieved by adding the quadratic term, 



5N ( A - 



(3.2) 



1 Note that in d = 3 this is only true if one considers 4x4 Dirac matrices. 



to Cs{ip, ip). We are then led to the interpolated model 



Cs = ipk (i $) - 5aipk^k - (1 - 5) i] ip k i/j k - -^-o- 2 + £ct,s , (3.3) 

where C c t,s is the part of Lagrangian density containing the necessary counterterms for renormalization, whose 
coefficients are allowed to be 5 and ?y dependent [13, [H| . As it is well known the 2 + 1-dimensional GN model is 
not renormalizable in the usual perturbative expansion, but is renormalizable in the 1/N expansion [25j ]. which has 
the property of modifying non-perturbatively the usual behavior under power counting. Though our renormalization 
procedure in the OPT expansion is more similar to a perturbative renormalization, this is not a real obstacle for our 
analysis. On general grounds it is always possible to calculate physical quantities in a non-renormalizable model, at the 
price of introducing new counterterms at successive orders, which simply means that the sensitivity to an implicit cutoff 
of the model is expected to be more pronounced than in a renormalizable theory [26] . Such a procedure is commonly 
and successfully applied in many effective theories, like e.g. typically in chiral perturbation theory (for a recent review 
of chiral perturbation theory with emphasize on the renormalization procedure see e.g. [27| and references therein). 
However, concretely in our case, at first order of the OPT expansion all the relevant quantities are actually finite 
(when using dimensional regularization as we do here), thus completely unambiguous. Next, at second order, that we 
also investigate to some extent in this paper, it turns out, that the only potentially non-renormalizable contributions 
to the effective potential actually vanish, such that only standard (i.e. mass, wave-function, etc) counterterms are 
necessary to cancel the divergences. It should be mentioned however that this is somehow an accident of using 
dimensional regularization, which therefore delays at most, i.e. to much higher orders, the necessary introduction of 
new counterterms in our case. A detailed account for renormalization of the GN3d model in the OPT up to second 
order, will be presented in Apps. B and C. 

From the Lagrangian density in the interpolated form, Eq. (|3.3p . we can immediately read the corresponding new 
Feynman rules in Minkowski space. Each Yukawa vertex carries a factor —iS while the (free) a propagator is now 
—iX/(NS). The LDE dressed fermion propagator is 



SAP) = p —- , (3-4) 



where 



77* = r? - (t) - cr c )5 . (3.5) 

Any quantity computed from the above rules, at some finite order in <5, is dependent on the parameter 77, which 
then must be fixed somehow. Here, as in most of the previous references on the OPT method, 77 is fixed by using 
the principle of minimal sensitivity (PMS). In the PMS procedure one requires that a physical quantity that is 
calculated perturbatively to some fc-th order in S, be evaluated at the point where it is less sensitive to this parameter. 
This criterion then translates into the variational relation 28] 



d$(fc) 
dr/ 



fj,S=l 



= . 



(3.6) 



The optimum value 77 that satisfies Eq. (|3.6p must be a function of the original parameters, including the couplings, 
thus generating non-perturbative results. 



A. The OPT Effective Potential 



The different contributions to the order-5 2 self-energy are displayed in Fig. [T] We can use these self-energy terms 
to evaluate the vacuum graphs contributing to the effective potential as shown in Fig. [2] 

In the sequel we make use of the following notations. The four-momentum P is given by P — (fb,p), where 
Pq = i(u> n — ifi), with Lu n — (2n + 1)tvT, n — 0, ±1,±2, are the Matsubara frequencies for fermions. The 
momentum integrals, when passing from Minkowski to Euclidean space-time, we here denote by 
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FIG. 1: Diagrams contributing to the self energy to order- J 2 . The thick continuous fermionic lines represent rj, dependent 
terms which must be further expanded, while the thin continuous lines represent r\ dependent fermionic propagators and the 
dashed lines represent the a propagator. Diagrams (a) and (b) (of order-<5 and order-5 2 , respectively) contribute with 1/JV, 
while diagrams (c) and (d) (both of order-5 2 ) contribute with 1/N 2 . Within dimensional regularization only graphs (b) and 
(c) are divergent. Tadpole graphs are not shown as they do not contribute to the effective potential nor to the counterterms 
(in d — 2 + 1) at the perturbative order we restrict ourselves to. 
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FIG. 2: Diagrams contributing to V e s/N to order-5 2 . The thick continuous fermionic lines represent r)* dependent terms 
which must be further expanded while the thin continuous lines represent r\ dependent fermionic propagators and the dashed 
lines represent the a propagator. The first (order-5 ) contributes with 1/N°, the second and third (order-<5 and order-<5 2 
respectively) contribute with 1/N. The fourth and fifth (both of order-(5 2 ) contribute with 1/N 2 . The sixth and seventh 
represent contributions due to the mass and wave function renormalization counterterms respectively. The last graph represents 
the zero point energy subtraction term. 



where all space momentum integrals are performed using dimensional regularization, d = 3 — e. The renormalization 
procedure (which is only necessary at order-(5 2 since all relevant quantities are explicitly finite at order-<5) is carried 
out in the modified minimal subtraction scheme (MS). 

The order-5 OPT effective potential is obtained from the first two diagrams shown in Fig. [3 and, using the previous 
Feynman rules for the GN in the OPT, it is given by 



V eS , s i (a c , 77) a 2 . 

■ = 0— ^ + 1 

N 2A 



(T) 



trln(,P - if) + Si 



(T) 



tr- 



V - o c 



P-r) + ie 



N 



(3.7) 



where AV}^ /N brings the first 1/N correction to the effective potential shown by the second diagram in Fig. [5] This 
is given by [2!| 



(a) 



N 



'2N 



f(T) 


\ S a (j7) I 


J>-r) + ie_ 


/ ^ 



(3.8) 



where the trace is over Dirac's matrices only 2 while the term S a represents the first contribution shown in Fig. [T]to 
the fermion self-energy, 



A 



1 



Iq K V + ie 

After taking the traces in Eq. (|3.7[) and rearranging the terms one obtains 



(3.9) 



N 



6 2X 



2 1 



(T) 



In (P 2 



8Ai 



(T) 



P 2 



2 The factor — 1 corresponding to a closed fermionic has already been taken into account |29ll . 
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A 2A 2 



(T) 



P 2 — ?7 2 + ze 



1 2 


V 




Jp 



-co Pn 



-1 2 



(3.10) 



Then, at finite temperature and chemical potential, one finds (see appendix A for the relevant integrals and Matsubara 
sums leading to this result) 

iV = 6 2\ + ^ + — h{a > b) + - l2{a > b) 

_ s vtv_^± M + rj3(0j + ^.^^ y + T j s(0j 5)] 2 + 5-^- [h(a, b)f , (3.11) 



2(2tt) 2 A^ 



2(2tt) 2 7V 1 



where we have defined the functions 



7i(a,6)=Li 2 [- e -( a - b )]+Li 2 [- 



-(o+6)l 



I 2 (a,b) =Li 3 [-e-( a - 6 >] +Li 3 [-e-( a+b )] , 



7 3 (a,6) = In 



l + e -(»-6) 



In 



1 + e-( a+b > 



(3.12) 
(3.13) 
(3.14) 



J 4 (a, 6) = sgn(/x) 



a In 



1 + e 



a+b 



1 + e a - b 



Li 2 [-e Q +VLi2he 



a—bl 



with a = l^l/T and b = \[J*\/T. 

The T — > limit for each of the elements appearing in Eq. p. lip are: 



YxmT 2 h(a,b) = --(\fj,\-\r,\y8(\n\-\r,\) , 

|iniP 3 / 2 ( a ,6) = i(H-| M |) 3 0(| M |-H), 

HmTJ 3 (o,6) = (|A*|-|»?|)e(|M|-N), 

hm T 2 h{a,b) = isgn( M ) (?7 2 — (J 2 ) 0(H - |r?| 



where #(|/i| — is the step function. 



(3.15) 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



IV. OPTIMIZATION AND NUMERICAL RESULTS BEYOND LARGE- TV 

Before proceeding to the specific d = 3 case, considered in this work, let us apply the PMS to the most general 
order-<5 effective potential, which is given by Eq. (|3.10|) . This exercise will help the reader to visualize the way the 
OPT-PMS resums the perturbative series. Setting 6 — 1 and applying the PMS to Eq. (|3.10[) we obtain that 



(T) 



N 



(T) 



V 

Pa 



p2 _ v 2 



di] 



(T) 



P 2_ rj 2 



p2 _ jj2 _j_ j £ / ^ 



(T) 



Po 



p2_ v 2 



(4.1) 



As one can see in App. A, Eq. (|A3[) . the last term of the above equation only survives when \i ^ 0. In the case /i = 0, 
Eq. (|4.ip factorizes in a nice way which allows us to understand the way the OPT-PMS procedure resums the series 
producing non-perturbative results. With this aim one can easily check that (at 5=1) 
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Then, when /x = 0, the PMS equation factorizes to 



(T) 



[77 - ct c - £ a (?7, T,fi = 0)] ( 1 + r\ 



dr/ 



p2 _ ^2 _|_ j £ 



(T) 



(4.2) 



P 2 — rj 2 + ie 



= 



(4.3) 



leading to the self-consistent relation 



rj = cr c + Yi a (r),T,n = 0) 



(4.4) 



which is valid for any temperature and number of space-time dimensions. In this way the OPT fermionic loops get 
contributions containing a c as well as rainbow (exchange) type of self-energy terms, like the first graph of figure [T] 
Note that when N — > oo, fj = a c and the large N result is exactly reproduced [HI. The mathematical possibility 



(T) 



p2 _ fj2 



(4.5) 



corresponds to the unphysical, coupling independent, solution discussed in Ref. fl3| . Note that in the d — 3 case one 
obtains 



£ a (?7,T, n) 



\ V \+TI 3 (a,b)} 



(4.6) 



For numerical purposes, taking A — > —A and defining A = 7r/|A|, one can consider the dimensionless PMS equation: 



71 ~ ° c + w iM + Th{a ' h)) \ i 1 + {lvl + Th{a ' h)) + W Ii{a ' b) ^ h{a ' h) 



= . 



(4.7) 



where rj, a c , T, and /i are in units of A. 



A. The T = and n = case 



Let us start by analyzing each of the different possible cases involving the temperature and chemical potential 
corrections. For T — and fi — we have that, to order-<5, 



N 



2A 3tt 



2(27r) 2 A^ 



(4.8) 



Note that (rj, a c ) — V^ s {—rj, —cr c ) and by the virtue of this symmetry we shall look for fj(a c ) only for a c > 0, since 
for a c < it is obvious that fj((T c ) — — fj{— a c ). 

Then, using A — > — A and A = 7r/|A|, one can write the free energy, at T — and /i = 0, as 



N 2tt + 3tt tt SttNA ' 1 Jj 

where the notation is consistent with the fact that we are only interested in a c > (in this case only rj > can recover 
the large- TV result as N — > oo as can be seen from the PMS solution, Eq. (|4.4p ). Then dV c s/da c = at a c — a c gives 
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a c = rj 2 /A . 

At the same time the PMS equation dV c « jdr\ — at 77 = fj gives the relation: 



(4.10) 



from where we then obtain the expression 



with the function F(N) defined as 



4NA ' 



(4.11) 



(4.12) 



(4.13) 



The above results then lead to the optimized value for the (dynamically generated) vacuum expectation value for the 
scalar field , also shown in [l8l |. 



A 



This result is contrasted with the large- N result in Fig. [31 



(4.14) 




A 



FIG. 3: The dimensionless minimum <r c (in units of A) as a function of A for T = \i = 0. The dashed line represents the 
N — > 00 result, while the continuous lines were produced by the OPT-PMS at order-5. The numbers beside the curves identify 
the value of N for each case. 

It is instructive to recall a similar result for the d = 2 where the large- N result is ct^(A) = M cxp(— n/X) 25], 
while the OPT result is rrf = ~ 1/(2N)], where A* = A[ l -1/ (2N)] [13]. The same happens here except 

that we have a factor 4 inside the function dependent on N, Eq. (|4.13[) . instead of a factor 2 found in the d = 2 
case. This is because of the 4x4 Dirac matrices considered in the three dimensional problem here. Thus, we have 
<j 5 c — <7^(A*)/[1 — 1/(42V)], recalling that ct^(A) = A = 7r/|A|. Finally, when we evaluate the thermodynamical 
potential in the sequel, it will be useful to consider the optimized free energy at its minimum, a c = a c , which is given 
by 



N 



A 3 



(4.15) 
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B. The TV and fj, = case 

Next, let us consider the case T ^ and (i = when the free energy can be written, using again A — > — A, and 
A = 7r/|A| as: 



N 



= -5A- 



.-■2- + - 

27T 37T 7T 



^T 2 Li 2 (-e-"/ T ) +T 3 Li 3 (-e-"/ T ) 



-(r7-cr c ) f) + 2Tb 1 + e 



-77/T 



?y + 2Tln 1 + e 



(4.16) 



where we have considered again the case a c > and 7/ > 0. From the result given in Eq. (|4.4p . 77 = cr c +S (rj, T, /j, = 0), 
we immediately obtain the self-consistent temperature dependent relation 



4NA 



+ 2Tln(l + e- f?/T ) 



(4.17) 



It is a simple matter to apply dV c s / da c — at a c — a c to Eq. (|4.16p to obtain 



V 

ac = A 



?7 + 2Tln(l 



-r,/T 



(4.18) 



which can be used in Eq. (|4.17j) to yield rj — aJ-(N) that, when inserted into Eq. ()4.18|) . allows to study the thermal 
behavior of order parameter (a c (T)) via the extremum of (77, a c , T), given by [18| 



A 



2T]n[l + e-^ T) ^ N)/T ] 



(4.19) 



From the above equation one retrieves the result <7 C (0) = AJ r (N)^ 2 , as obtained in the previous subsection. The 
critical temperature T c for chiral symmetry restoration is obtained by requiring that a c (T = T c ) = 0, which gives the 
result 



rpo 



A 



2\n2J : -(N) 



(4.20) 



This analytical result is shown in the Fig. [4j Note that a numerical application of the PMS to the OPT effective 
potential, using Eqs. (|4.7|) and (|4.16|) . exactly reproduces the analytical result Eq. (|4.20j) . Recall the similarity with 
the d = 2 case [l3|, where the LDE result was obtained from the large N result by the replacement: A — > A* = 
A[l — l/(2N)], just like the same as in obtaining the scalar field vev, as discussed in the previous subsection. 

Note also that, contrary to the d = 2 case, our prediction for T c is always greater (for finite N) than the large-Af 
prediction. The transition is found to be of the second kind, as illustrated by Fig. [5] 

C. The T = and fj, ^ case 

Let us now consider the case T = and fj, 7^ 0. From the general expression of the effective potential at the first 
OPT order, Eq. (|3.10p . and using the T — > results shown in Eqs. (|3.16p - (|3.19p . we find that the chemical potential 
dependent effective potential is given by 
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FIG. 4: The OPT critical temperature (at /j, — and in units of A) as a function of A. The dashed line represents the N — > oo 
result, while the continuous lines were produced by the OPT-PMS procedure at order-5. The numbers beside the curves identify 
the value of N for each case. 




FIG. 5: The large- A'' (dashed line) and the OPT (continuous line) predictions for o c (T) at N — 3. All quantities are in units 
of A. The figure displays a continuous, second order, transition line. 

For (i = 0, obviously, Eq. (|4.2ip reduces to Eq. (|4~8|) . 

To evaluate the critical value /i c for chiral symmetry restoration, it is sufficient to compare the values of the effective 
potential at the minimum V c g(a- C , T = 0) with its value for fi ^ for a c = 0. This is from the same line of reasoning 
employed in the d = 2 case discussed in Ref. [l3j]. In this case, we obtain the point where the two minima of the 
effective potential, at a c — and at a c ^ and \i = /i c coincide, i.e., there is a value /i c which satisfies 



V& [a = 0, M = Mc , T = 0) = Kff (*c M = 0, T = 0) . (4.22) 
It is a simple algebraic exercise to calculate both members of this equality. We first obtain from Eq. J4.21j) that 
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Then, to evaluate the right hand side of Eq. (|4.22[) we use Eq. (|4.8[) together with the relation between a c and 77 in 
Eq. P~T2"|) , i.e., ct c = 7?/[J"(iV)], which gives 



where the last simplification arises from using the definition of F(N), Eq. (|4.13p . and noting that rj = —tt/[XT(N)]. 
Note thus that V^g (ct c , h = T = 0) has formally the same simple expression as the leading order one, except of course 
that it includes non-trivial 1/N corrections via the explicit expression of fj. Now we may compare Eqs. (|4.23|) and 
(gHI to finally extract fi c [l|], 

= (^4Tr¥Y / \ (4-25) 



f(N) \ 16N A 
where we used again A = vr/ 1 A| . For AT = 3, we find from Eq. (|4.25[) the solution 



^ ~ 1.06767 , (4.26) 



which agrees with the numerical results obtained in the next subsection. 



D. The T/0 and ft ^ case 



Finally, turning now for the case of both finite temperature and finite chemical potential, we obtain the full phase 
diagram of the three dimensional GN model. The numerical application of the PMS shows how the phase diagram is 
qualitatively and quantitatively affected by finite A" corrections. Figure [6] shows the situation for A" = 3. The CSB 
region is augmented with respect to the large- A^ predictions, when expressed in units of our reference scale A = 7r/|A|. 
This is clear also from our results for a c , T c and /i c . This appears at first sight in contrast with the lattice results 
which show a decreasing of the CSB region at finite N, as compared to the large- N results. However, whithin 
our approximation, the increase of the size of the CSB region is rather small, being about 5% for N = 3, while the 
increase would be of only about 2% for the N = 12 case considered in Ref. which in turn predicts a decrease of 
about 10%. One should moreover note that in Ref. [13] the authors present their results for the phase diagram with 
quantitites normalized by the scalar vacuum expectation value obtained from the lattice simulations. In our case it 
means that from Eqs. (|3~T4"|) . (|4~20f and (j4~25|) . taking N = 3 for instance, T c /a c ~ 0.661 and [i c ja c ^ 0.897, while 
for N = 12, T c /a c ~ 0.707 and fJ. c /<7 c ~ 0.976, which approximately agrees with the results presented in Ref. [l2j 
within their level of precision. Note also that our results agree reasonably with recent analysis of the 2 + 1 GN model 
from exact renormalization group methods (24| (at least within the errors quoted there, and for the case of vanishing 
wave function of the auxiliary scalar field, which is the appropriate comparison to our analysis). Thus the reduction 
or increase with respect to large N results appears to be just a matter of scaling. In the present work we consider the 
scale defined as A (the vacuum expectation value of the scalar field at large- N) as more appropriate to present the 
results. In a previous work [l3| . we have considered the same model but in 1+1 dimensions using exactly the same 
approximation. There, our result for the phase diagram showed a drastic change concerning the size of the CSB which 
was about 30 % smaller than the one produced by the large- A^ approximation for A" = 3. Now, in 2+1 dimensions, it 
turns out that the CSB region predicted by the OPT is very close to the one predicted by the large- N approximation. 
In summary, in the fixed normalization scale used (A) , it looks like the OPT predicts a drastic decrease in the size of 
the CSB region in 1+1 dimensions whereas in 2+1 dimensions it seems to support, at least at lowest OPT order, the 
CSB size (as well as the numerical values of a c , T c and // c ) predicted at large- N. We shall see in next section that 
a partial investigation of higher OPT order corrections (for T = /j, = 0) indicates a good stability of these first order 
results. On the other hand, the nature of the transition line predicted at large- A^ in 1+1 dimensions is unaffected by 
the OPT at order-(5 while it drastically changes in 2+1 dimensions, as we now start to discuss. 
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FIG. 6: The large-TV (dashed line) and OPT (continuous line) predictions for the phase diagram at TV = 3. All quantities are 
in units of A. The black dot indicates the position of a tricritical point, located at T tcI ~ 0.251 and im, ci ~ 1.029. Below this 
point the transition is of the first kind while above the point it is of the second kind. 

Concerning the nature of the transition lines shown in Fig. recall that the large- N approximation predicts that 
it is of the second kind everywhere, except at T = where it suddenly becomes of the first kind at a critical value 
of /i = fj,^ = A. Using lattice Monte Carlo simulations for the GN model, Kogut and Strouthos [ll| concluded that, 
for N = 4, there should be a tricritical point on the section of the phase boundary defined by T/T^ < 0.230 and 
n/Hc — 0.970. Our evaluations predict that the observed first order transition at T = spreads out through the 
transition line until it reaches a tricritical point at T tcr ~ 0.251 A and /iter — 1.029 A (for N = 3). 

Figure shows an envelope of curves that displays how the abrupt first order transition at T = becomes smoother 
as the temperature increases. One observes that the discontinuity gap becomes smaller as the temperature approaches 
the tricritical value, T tcr = 0.251 A. 

It is also important to analyze the occurrence of metastability lines related to the first order transition line (T < T tcr ) 
in Fig. [51 With this aim we offer Fig. [5J where the dashed line joining the tricritical point, P t , to point A corresponds 
to the appearance of a minimum at a c = 0, whereas the continuous line joining points Pt and /i c refers to the first 
order transition line. Finally, the dot-dashed line joining points Pt and B refers to the vanishing of the minima that 
occurs away from the origin. It is important to note how small the metastable region A — P t — B is. As a matter of 
fact, point A occurs at a value of /i which is about 2% smaller than /i^, whereas B occurs at a value which is about 
3% greater than \x c . In d = 2 these values are of about 30% (see Ref. @). 

The reader may visualize the three situations shown in Fig. [5] by examining the form of the free energy shown in 



We are now in position to perform a more physical interpretation of our results by examining other relevant 
thermo dynamical quantities. The thcrmodynamical potential, f2(/x, T), for instance, is related to the free energy at 
its minimum, a c — d c . It is given by 



Note that we have defined this quantity in terms of the optimized free energy. This is an important remark because 
we are considering different physical quantities and one could wonder which one to optimize. Here, our choice is to 
optimize the free energy since all other thermodynamical quantities may be obtained from it. It is usual to normalize 
the thermodynamical potential by subtracting a "bag" term, B, given by B = 0(0, 0), so that the pressure as well as 
the energy density vanish at T = and /j, — 0. In view of Eq. (|4.15|) the bag term is simply given by 
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V. THE LIQUID-GAS PHASE 



tt(n,T) = Ves(Tj,ff c ,H,T) . 



(5.1) 
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FIG. 7: The order parameter, a c , as a function of /x for different temperatures. The continuous lines represent the OPT results 
for N = 3 and the labels on the figure represent the different temperatures: T a = 0.050, Tt = 0.100, T c = 0.150, Td = 0.200 and 
T e = 0.250. All these values are smaller than the tricritical value T tcr = 0.251 and the associated curves clearly display first 
order transitions. All quantities are given in units of A. For reference, the figure also shows the large-iV result (dashed line) 
for T = 0.150. 
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FIG. 8: Part of the phase diagram that corresponds to the metastable region (reproduced here from [Hj]). The dashed line 
joining points Pt and A refers to the development of a minimum at the origin. The continuous line linking the tricritical point 
to C (T = and /i = fj, a ) is the first order transition line, while the dot-dashed line joining Pt to B is the second metastability 
line and signals that the minima that occur away from the origin have disappeared. 



B = -^Nf (5 - 2) 

Then, the normalized thermodynamical potential is just fijv(/i, T) = T) — B. At the same time, the (normalized) 
pressure is given by P(fi,T) — —Q,n([i,T) from which one may obtain the density 
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FIG. 9: The shape of the free energy corresponding to the metastable region. The left panel shows the situation corresponding 
to the dashed Pt — A line shown in Fig. [8j while the middle and right panels correspond to lines P t — C and P t — B respectively. 



dP 

(5.3) 

and the entropy density 

dP 

8- w . (5.4, 

Finally, the (normalized) energy density is given by £ = —P + TS + pp. Recalling that, due to the gap equation, 
dP/da c = and that, due to the PMS equation, dP/dfj = one obtains the density 

f)T\ T 3 , (fj-a c )fj rrT AryS XT 4 

P = -—/i,, - ~h,» + -^—Th,» {r, + Th) h„ - j^hh,, , (5.5) 

where Jj i(L1 = dli/dp, and Ii, i = 1, ... ,4, are given by the function in Eqs. ([312)) - (|3T5| . At the same time one 
obtains that the entropy density is given by 

o VT T fjT 2 T 2 T 3 , (v-<tc)fi T , (fj-a c )fj mT 

O — —Z ii 1\ t — 6 12 H,T H -<3 H 1 J3,T 

TT TT TT TT IT TT 

\f) 2 XT 3 
- WfN ^ + Th) {h + TI ^ T) - (2^N^ + Th ^ ' (5 - 6) 

where I^t = dli/dT. 

Having all the above quantities, we can now analyze, for instance, the phase diagram in the physically more 
accessible P — 1/ p plane as shown by Fig. [TDJ This figure displays one of our most important results which indicates 
that a mixed "liquid-gas" phase, previously unknown to exist within this model, develops at low pressure values. 
Three isotherms are shown and the one corresponding to T — defines the edge of the region accessible to the system. 
The dotted line above the tricritical point is just the mapping of the corresponding second order transition line in the 
T — p, plane. Note that in the P — 1/ p plane the first order transition line displayed in the T — p plane (corresponding 
to T < Iter) splits in two parts corresponding to the value of the density at the two degenerate minima which produce 
identical pressure values. Not being able to determine the existence of the mixed phase, Kogut and Strouthos argued 
that the liquid-gas transition was either (i) extremely weak, (ii) very close to the chiral transition, or (iii) not realized 
in this model [11(. Our result shown in Fig. [TOl suggests that their second hypothesis was the correct one. Moreover, 
these authors had used N — 4, that was the smallest number allowed by the hybrid Monte-Carlo algorithm used in 
their simulations. With this number the tricritical point appears at an even lower value of P (actually, for N — > oo 
it happens at P — 0) and was consequently harder to be detected within their approximation. 

Figure [14] shows a detailed view of the "liquid-gas" phase seen in Fig. [TO] displaying how an isotherm whose 
temperature is smaller than the tricritical value crosses the mixed phase region. The horizontal line in the "coexistence" 
region was drawn by connecting the value of the pressure at the boundaries, corresponding to "mixed" states. This 
picture observes the Maxwell construction which derives from the equality of the chemical potentials at the edge of 
the two phases. 

Figure [TOl shows £/T 3 and P/T 3 as functions of the temperature for p = 0. Note that for high temperatures 
£/T 3 -> -3C(3)/tt ~ 1.14, while P/T 3 -> -3C(3)/(2tt) ~ 0.57 (where C(3) ^ 1.202), as one can guess by looking at 
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FIG. 10: The phase diagram in the P-l/p plane for iV = 3. The thick continuous line is the T — isotherm which limits the 
region accessible to the system. The region labeled by N/A is not accessible. The dotted line is the mapping of the second order 
transition line. The chiral symmetric region (CSR) corresponds to the "liquid" phase while the region where chiral symmetry 
is broken (CSB) corresponds to the "gas" phase. The LG region starting at the the tricritical point, Pt, is limited by first 
order transition lines and corresponds to the mixed "liquid-gas" phase. Point C (P = and p ~ 0.197 A 2 ) corresponds to 
T = 0,/i = He- The isotherm represented by the dashed line corresponds to the tricritical temperature, T tcr = 0.251 A while 
the dot-dashed line represents the isotherm corresponding to T = 0.397 A. The pressure, P, is in units of A 3 while the density, 
p, is given in units of A 2 . 

the equations for £ and P at /z = 0. In those high temperature regimes (T > T c ), we have a c = and fj — > 0. At high 
temperatures both curves are symmetrical with respect to the numerical value ~ 0.85. 

Let us now check for the presence of latent heat, which is inherent to first order phase transitions. We do this by 
examining the energy density as a function of the temperature. For this, one chooses a value of \i that corresponds 
to the first order transition, like, for example, any \± such that fj, c > p, > /iter- Recall that for the case N = 3, 
/iter = 1.029A and \i c = 1.067 A, so we choose, without loss of generality, the value \i = 1.040A. One then expects 
to see a discontinuity in the line corresponding to £{T) at T = T c (/i = 0.140A) = 0.194A. This is indeed the case as 
shown in Fig. 1131 The same figure shows the large- N result, where the discontinuity happens only at T = 0, which 
can be understood by recalling that within this approximation the first order phase transition happens only at the 
point T = and n = 1.000A. 



Let us now investigate the order-5 2 contributions that are given by the three-loop graphs shown in Fig. O Actually, 
a complete evaluation of these graphs at finite T and fi turns out to be very cumbersome, so we shall restrict ourselves 
in the present work to the T — and [i = case which is more tractable. This will at least allow us to have a 
reasonable quantitative estimate of the expected higher order corrections to our previous order-J results. We plan to 
tackle the calculation of the full T and fj, dependence in a future work. 



VI. ORDERS 2 RESULTS AT T = AND p = 



A. The order 5 2 three-loop contribution to the free energy 



The total order S 2 three-loop contribution can be easily extracted from Ref. [2!| and reads 




(6.1) 
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FIG. 11: Detail of the liquid-gas phase in the P-l/p plane for iV = 3. The isotherm represented by the dashed line corresponds 
to a temperature T = 0.194 A that is smaller than the tricritical temperature, T tcr = 0.251 A . The two dots are joined by a 
straight line in a Maxwell construction. The pressure, P, is in units of A 3 while the density, p, is given in units of A 2 . 




where T,i(p,rf), i — b,c,d, correspond to panels (b),(c), and (d) (second, third and forth diagrams, respectively) in 
Fig.ffl 

The most complicated contributions arise from the first and second terms of the above equation (or equivalently 
corresponding to the third and fourth graphs of Fig. [2]) since in this case the self energies depend on the momentum 
p. After taking the traces, etc and using dimensional regularization with d = 3 — e, the corresponding integral to be 
evaluated reads 
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FIG. 13: The energy density, £, as a function of the temperature for /i = 1.040 > /iter = 1.029 (for TV = 3). Both quantities 
are in units of A. The continuous line is the OPT result and shows the presence of latent heat signaled by the discontinuity at 
T c (p = 1.040A) = 0.194A. The large- TV result is represented by the dashed line and the discontinuity happens at T = 0. Both 
£ and T are in units of A. 
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Introducing next appropriate Feynman parameters to disentangle the different momenta integrations, we obtain the 
following expression 
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where M is the arbitrary MS renormalization scale, 



?(7) = [7(l-7)r 3/2+£/2 , 



(6.4) 



and 
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(6.5) 



The evaluation of the final integrals over Feynman parameters in Eq. (|6.3[) is rather technical and details of this 
calculation are left to appendix B. One arrives at the final result for the third and fourth three-loop diagrams of Fig. 
[2] as given by 
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where X ~ 1.63669 is a numerical constant obtained from the integrations. Next, the self-energy entering the last 
term of Eq. (|6.ip . corresponding to the last diagram shown in Fig. [I] that leads to the contribution to V e g shown by 
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the fifth diagram in Fig. [2j it is a tadpole like graph that can be easily evaluated with the standard Feynman rules. 
It is finite and gives the contribution, 
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At the three-loop order, the free energy contains divergent terms that lead to the 1/e term shown in Eq. (|6.6p . 
The renormalization is performed as usual, by introducing the appropriate counterterms (mass, wave function, etc). 
Actually the perturbative two-loop fermion self-energy in Fig. 14 exhibits divergent terms of non-renormalizable 
kind, with higher power of momentum dependence, as expected since the model is not perturbatively renormalizable. 
Although it would not be a problem in principle to treat those new divergences with appropriate counterterms, 
similarly to what is done in other effective theories, it turns out that these non-renormalizable counterterms do not 
contribute to the three-loop free energy in dimensional regularization, so that only standard mass, wave-function and 
vacuum counterterms are needed in practice to render the effective potential at 0(S 2 ) finite. Thus the renormalization 
as performed in the MS scheme after dimensional regularization does not introduce new parameters at the three-loop 
order, for the quantities we are interested in. The arbitrariness of the final physical result for the three-loop effective 
potential simply takes the standard form of a (logarithmic) dependence on an arbitrary renormalization scale. The 
details of the calculation of these counterterms contributions are discussed in App. C. Now, by adding all contributions, 
including the finite ones arising from the counterterms (see App. C), the renormalized three-loop effective potential 
can be cast into the form: 
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B. Optimization results at T = and — at order 8 2 



From the result for the 0(X 2 S 2 ) three-loop contribution, Eq. (|6.8[) , we can now perform the S expansion and PMS 
optimization at this next order, limited however here to the special case T = and fi — 0. Following the same 
line of reasoning as in Sec. IV. A above, care is to be taken by noticing that other S 2 terms are generated by the 
appropriate expansion of 77*, as defined in Eq. (|3.5j) . within the first order 0{8\) terms. Explicitly, after having 
redefined A — > — tt/A as previously, we arrive at the expression, 
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(6.9) 



The remaining is simply an algebraic exercise to apply the PMS procedure to Eq. (|6.9p . Similarly to the first order 
case, the gap equation, dV c s/da c — at a c — a c defines a c as function of fj, while the PMS equation dV c s/dr] = at 
77 = fj gives a further relation between a c and fj. However, at this next order, both relations are more complicate, in 
particular the PMS equation is non-linear and involves a ln(?7) term. It is most convenient to use the gap equation, 
which gives 
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generalizing Eq. (|4.10[) . at second 5-order, into the PMS equation, which defines an equation depending only on fj(N) 
(note that the dependence upon A can be simply factored out, e.g., by rescaling fj in units of A). For instance, for 
the particular case N = 3, which is just sufficient for our illustration, one obtains the PMS equation as 
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where we have defined, for convenience, the dimensionless mass parameter fj = r)/A. Eq. (|6 . 1 1 [) can be solved 
numerically, to find a non-trivial value of fj. For this we have to set the arbitrary renormalization scale M in (|6.11[) . 
originating from the logarithmic dependence in Eq. (|6.9|) . to some appropriate value. A physically natural choice 
is to set M = A, that corresponds to the basic scale and scalar vacuum expectation value in the large- N limit. An 
interesting feature of the optimization result in the present case, is that the presence of the In (77) dependence in Eq. 
(|6. 1 1|) largely reduces the number of optimized solutions (usually a drawback of the PMS). Indeed, one sees in Eq. 
(|6.11[) that, without the logarithmic term, four different non-zero (real or complex) 77 solutions would occur 3 . In 
contrast, we find here numerically a unique (real) solution (for M = A and N = 3): 

77 - 0.867 A. (6.12) 
Next, we can just plug in this result into Eq. (|6.10[) , to obtain 

of ~ 1.1719 A , (6.13) 

which appears to be very close to the first order result, of ~ 1.1901A, obtained from Eq. (|4.14[) for N = 3. One 
may however question if this is an artifact of our choice of arbitrary renormalization scale A = M. Though this value 
of the scale appears to be very natural, it is easy to study the impact of varying it in a reasonable range around 
this value, like is sensible to do in similar renormalization scale dependence studies in other theories (and which 
give a rough estimate of higher order corrections). We can, e.g., vary M in the range 0.50A < M < 1.50A, which 
correspondingly changes o c from values 1.16A < o c < 1.38A (note that for M values too far apart from M — A, 
Eq. (|6.1ip has no real solutions). From this we conclude that the higher 5 2 order corrections for T, /1 = to the 
previous analysis are quite small, though their scale dependence is apparent and can become non negligible. Actually, 
it is clear that the scale dependence is essentially determined by the relative size of the hi(rj/M) term with respect to 
constant terms in the optimization Eq. (|6.11[) . which turns out to be moderate 4 . On general grounds, renormalization 
scale dependence is expected to be rather pronounced at lowest loop orders and damped by higher order perturbative 
contributions [26l |. But in our framework the three- loop expression, Eq. (|6.8p . is the very first perturbative order 
at which renormalization, and thus scale dependence, occur for the effective potential, so that the behavior is more 
similar to a one-loop quantity with respect to renormalization scale dependence. For the natural scale choice M = A, 
however, these second order corrections for cf c value in Eq. (|6.13p are very small, about only 1.5% of the first order 
value, so in other words the OPT expansion seems to converge rather quickly. 

VII. CONCLUSIONS 

In the present paper we have applied the alternative analytical non-perturbative optimized perturbation theory 
(OPT) approach, through which one can easily include finite N effects, to a four-fermion theory with discrete chiral 
symmetry represented by the massless Gross-Neveu model in 2+1 dimensions. We have then evaluated the free energy, 
or effective potential, for the model at both finite temperature and finite chemical potential. The evaluation of the 
optimal value for this quantity has allowed us to derive and study in detail other thermodynamical quantities, such 
as the pressure, density, entropy and energy density. 

Our main results in this paper include the analytical derivation of expressions, going beyond the standard large- N 
results, for the scalar field vacuum expectation value, the critical temperature, and for critical chemical potential 
related to chiral symmetry breaking/restoration. We have also demonstrated the existence of a tricritical point, not 
seen in the large- N approximation and the corresponding existence of a mixed chiral restored-broken phase in the 
system for finite values of chemical potential and temperature. Concerning the phase diagram we recall that in a 
lattice simulation, Kogut and Strouthos [llj have predicted the existence of tricritical point in the T — /_t plane that 
is missed by the large- N approximation. However, within the numerical precision of their simulations, those authors 



3 One may wonder if the peculiar choice of the arbitrary scale: M = fj, thus canceling the logarithmic term, would not re-introduce the 
PMS multi-solution problem. But this gives no consistent optimal solutions since all other fj solutions contradict this value of fj. This 
incidentally shows that Eq. (16.111 ) does not always have real solutions for any values of M. 

4 It is interesting to note also that from Eq. (16. 10 ^ o- c (rj) has a minimum, a c ~ 1.16 A, at 77 ~ 0.94A, independent of the scale M, thus 
a c cannot reach values below ~ 1.16A: when varying the scale M to values lower than A, the optimal solution fj of Eq. J 6- 1 1 ft decreases, 
so that iic(fj) will start to increase again. Also, the unwelcome singularity of <r c for r\ = A/2, apparent in Eq. Il6.10j l. is not a solution 
of Eq. 16.111 . so this singular value is never reached by the consistent PMS optimal solution fj. 
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were unable to give its exact location. Here, we have not only confirmed the existence of such a point but have also 
demonstrated how it can be located for any value of N. The formalism employed in this work has allowed us to easily 
draw the first order transition line together with the metastability lines. This exercise showed that the metastable 
region is rather small which possibly explains why its observation was not possible in Refs. (ill . H3 |. Going to the 
P — l/p plane has generated another important result for the GN3d model. Namely, the prediction of a "liquid-gas" 
phase transition, so far unknown to exist in this model. All these results drastically change the large-A^ picture of 
the GN3d phase diagram in which only a "superconducting" phase and a "normal" phase appear [J]. Although in 
our study the complete T and /i dependence could only be studied at the first non-trivial order of the S expansion, 
we were able to estimate the higher order-5 2 corrections at least at T = and /j, = 0. These corrections turn out 
to be reasonably small (about only 1.5% for a natural choice of renormalization scale). This also indicates that all 
others OPT order-<5 2 results are not expected to distance too much from the computed order-<5 results obtained for 
the GN3d model. So we may argue that all our phase transition and thermodynamical results obtained in Sees. V 
and VI are likely to remain qualitatively unaltered by higher order corrections in this framework. 

We notice that some of our quantitative results (most notably those displayed by Figs. 3 to 6) appear at first sight 
different from the overall behavior obtained with some Monte Carlo simulations [12| that predict a reduced size for 
the CSB region, while our results display an increase of the CSB region with respect to the large N (equivalently 
mean field) results. As discussed previously, as far as we can compare this seems to be only an artefact of the different 
reference scales chosen to express the critical quantitites (e.g. the authors of Ref. [12] express the phase diagram in 
terms of the scalar field vacuum expectation value, while we express all the relevant quantities in units of the fixed 
reference scale (A). But we note also that even for our choice of reference scale the increase with respect the large- N 
results is rather very small. This is opposite to what is seen when the same approximation is applied to the 1+1 
dimensional GN case [l3| where a large decrease of the CSB region is observed compared to the large- TV (mean-field) 
results. In this respect, the OPT applied to the 2+1 dimensions GN model leads to results closer to the mean-field 
approximation, though the character of CSB/CSR transition line is qualitatively very different. 

As usual, within the OPT, all the large- N results are exactly recovered from our expressions by taking the limit 
N — » oo. In addition to the explicit thermodynamical analysis presented here, we have also shown in details, in 
Apps. B and C, the renormalization of the effective potential up to 0(6 2 ) in the OPT. A comparison to the case of 
renormalization within the 1/N expansion is also briefly provided. In the case of the main expressions derived in this 
work, all results to 0(6) are shown to be finite within the dimensional regularization in the OPT. Divergences start 
to appear at second order. At this order the effective potential can be rendered completely finite just with standard 
counterterms. Possible extensions of the present work are presented in Ref. [l 81 ] . 
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APPENDIX A: SUMMING MATSUBARA FREQUENCIES AND RELATED FORMULAS 



In this appendix we give the results for the main integrals and Matsubara sums appearing along the text. Following 
the general procedure for evaluating these sums [l9j and using dimensional regularization in the MS scheme, the 
momentum space integrals are written as 

f d 2 p fe iE M 2 Y /2 r d^p 

where d = 3 — e, M is an arbitrary mass scale and — 0.5772 is the Euler-Mascheroni constant. 
For the general case of T ^ and \i + 1 0, we then obtain, for example, that 



/ (T) ln(P 2 -, ? 2 + i£ ) = ^ + ^T 2 |Li 2 



_ e -(M-M)/r 



Li 2 



„-(M+H)/r 



T 3 r. 
+ — {Li 3 



_ e -(M-M)/T 



Li 3 



_ e -(\v\+M)/T 



(Al) 
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(A3) 



where sgn(/z) is the sign function and Li 1/ (z) is the polylogarithm function and it is defined (for v > 0) as [30| 

fe=i 

APPENDIX B: EVALUATION OF THREE-LOOP FREE ENERGY GRAPHS AT T = AND p = 0. 



In this appendix we give some technical details on the evaluation of the two complicated contributions to the 
three-loop free energy. These are the third and fourth graphs of Fig. 2 respectively, which give the expression Eq. 
(|6.3[) . There are probably different ways to perform this integral. For convenience we choose to first integrate on the 
Feynman parameter 7, which can be done analytically, with the result: 
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where 2-F1 is the hypergeometric function, with z = a(l — a)/[/3(l — /3)]. Next we can use some well-known properties of 
the hypergeometric function [3(|, relating 2 F\ [...; z] to 2-Fi[---; 1— z], and the power expansion in z of the hypergeometric 
function as 



iFi\p,q,r;z] = 



(p)k(q)k 
(r) k k\ 



(B2) 



where the (p)k etc are binomial coefficients. These manipulations allow to factorize explictly the a and [3 integrals in 
the simple form: 



^ ' T(2+p + q) 



(B3) 



with x = a, p. with various possible values of p, q (which in the framework of dimensional rcgularization is valid for 
any values of p, q, since it will eventually exhibits the location of the different poles in e = 0). Next it is just a matter 
of systematic algebra, using for instance Mathematica [3l|, to perform all these simple integrals, re-summing the 
resulting power series to all orders, while picking up the divergent and finite pieces we need. Actually the divergent 
part only occurs only in the first two terms of the power expansion in z of the 2 F± function, which is simple to extract 
analytically. In contrast, the finite part results from contributions of all terms, which can be extracted nuerically as 
the corresponding series converges fastly. Expanding all factors contained in Eq. (|Blj) for e — > finally yields 
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where X ~ 1.63669 is a constant coming from the integrations. The renormalization of the bare three-loop result, 
Eq. (|B4|) . is performed in the next appendix by introducing appropriate counterterms. 
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APPENDIX C: RENORMALIZATION OF THE THREE-LOOP FREE ENERGY 

Here, we give details on the renormalization procedure used. At order-5 2 new counterterms, which do not posses 
the form of the original tree-level Lagrangian, appear. This clearly is a manifestation, within our OPT framework, of 
perturbative non-renormalizability. These counterterms originate from higher momentum-dependent divergences of 
some of the two-loop fermion self-energies shown in Fig. [T] and then potentially enter as contributions to the three-loop 
effective potential (or cquivalently the free energy). Since the model is renormalizable in the 1/N expansion, it is clear 
that these non-renormalizable counterterms are perturbative artifacts, which can be seen indeed to disappear in the 
corresponding 1/N quantities, as we will briefly show below. But since we are not using the 1/N expansion explicitly 
in our framework we have to deal to some extent with these non-renormalizable terms. As we shall see below, however, 
the non-renormalizable contributions to the effective potential actually vanish, such that only standard mass, wave 
function, and zero point energy counterterms are necessary to cancel the divergences (though this feature is simply 
an accident of the effective potential calculated at the second perturbative order) . 



k+p-q 
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p-q 
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j t j „ \ k+p-q. 



q P P q \ k / P 

q-k 



massdiv WFdiv 

FIG. 14: Momentum dependent two loop graphs contributing to the fermionic self energy to order-<5 2 . £f,(p) is of order 1/N 
while S c (p) contributes with 1/N 2 . The other two graphs represent the mass and the wave function counterterms as indicated 
in the figure. 

At order-<5 2 we have three two-loop contributions to the fermionic self energy as shown if Fig. [TJ From those, only 
the two displayed in figure [14] are divergent and will be evaluated here in detail. First, note that by choosing the 
momenta routing in an appropriate manner as indicated in Fig. [TJ] one only needs to evaluate the first graph since 
the closed fermionic loop contributes with — AN. One may explicitly verify that + S c (p) = [1 — l/(4^V)]Sj(p). 

With the LDE Feynman rules at T — and /j, — one has 5 

;y ( s_ ( W-* A V [ til dd<1 jCtoL f W + f l) m+^-4)+V\ \ 

h[P> [ ' \5N ) J {2n) d (2TT) d {q 2 - V 2 ) {(k 2 ^^) [(k+p-q) 2 ' [ J 

Here all integrals are done in arbitrary dimension d = 3 — e. Then, after taking the trace we have 

y( ,_ , 2 4A 2 2 f d d k d d q {j+rf) [P + k(p-q) + r 1 2 } 
b[P) N J (2ir) d (2w) d (q 2 - rf) (k 2 - rf)[(k + p - q) 2 — rf] 1 ' 

Now one can introduce one Feynman parameter (a) to merge the two k dependent propagators. Then, one performs 
a Wick rotation (using dimensional regularization in the MS scheme) to carry on the integral over k. This yields 



4A 2 r 1 
E6(P) = - ^(4^3/2 (^M 2 Y' 2 T{e/2 - l/2)(2 - e) jf da[a(l 



«)] 1/2 - 



s/2 



5 The minus sign due to the fermionic loop has already been taken into account in Eq. IClt . 
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7 (27r) rf (q 2 E + V 2 )[(p-q) 2 E +^]- 1/2+e/2 ' 1 J 

where 7y 2 = 77 2 /[a(l — a)]. Now, one can use a second Feynman parameter (/3) to merge the two final propagators 
obtaining 



7V(4tt) 

where 



E 6 (p) = _^_*£L( e 7*M a ) e r(e - 1)(2 - e) f dadf3[a(l - a)] 1 ' 2 ^' 2 ^ - /3 )-3/2+e/2 W P)+v] ^ (C4) 



<^ = -/(l-/3)/3 + r? 2 /3 + »? 2 (l-/3). (C5) 

Note that in the above equations we have already returned to Minkowski space. Let us now examine the type 
of divergences and corresponding counterterms emerging from these two-loop fermion self-energy expressions. The 
counterterms appear in the Lagrangian density C s ct 



C 5 ct = {b k [i0 A s (rj)] ^ + B s {7 1 )^ k , (C6) 

which is to be added to the original Lagrangian density, Eq. (|3.1[) . The perturbative order at which the required 
counterterm contributions enter first can be readily found from Eq. (4.5): 

^jf-(w)=i J ^trln[ff(l+AV»7(l + B«] . (C7) 

where, implicitly, those counterterm coefficients are of order (5A) 2 (since at lowest XS order, all calculations are finite 
in dimensional rcgularization as we saw in section III and IV). These counterterms are extracted from the $ and 77 
dependent terms in Eq. (|C4[) . For the f term for instance: 
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(C8) 

where the only pole is contained in r(e — 1). One can then expand in e and perform the (finite) integrals is a and 
to obtain a counterterm as ffA s (r),p) where 

Similarly the other counterterm can be extracted from the r\ dependent term in Eq. (|C4[) : 



[£ b (p) + £ c (p)] massdiv = -^± ( 1_ 4^) (^ E M 2 rr[e-l](2-e) j\adP[a(l~a)]i-Hl-/3)-i+i 



(CIO) 

The same type of manipulations that lead to the pole in Eq. (jC9[) above fixes the other counterterm as r]B s (r],p) 
where 

s, <«) = j2 (4M 1 -^)s(- ,,2+9 *' 2 )- < C11) 

Note in Eqs. (IC9|) and (|C11[) the extra p 2 dependence which is the manifestation of the perturbative non- 
renormalizability. 
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FIG. 15: Next-to-leading order in the 1/N expansion graphs contributing to the lermionic sell energy. 

As mentioned above, all such non-renormalizable terms should be absent in the renormalizable 1/N expansion. 
More precisely, one can check this at the next-to-lcading 1 /N order where the equivalent contribution to the fcrmion 
self-energy is shown in Fig. [15] (which corresponds to an infinite sum of perturbative one-loop graphs). We can easily 
obtain the structure of the divergent terms, replacing expression (|C9[) and (|C11[) . in dimensional regularization as: 

^KH ■ (ci2) 

where we used expression 2.12b of Ref. [2f| for the d — 2 + 1 resummed <r-propagator. The expression (|C12[) 
is consistent with the calculation of wave-function and mass counterterm as performed in Ref. [25j with cutoff 
regularization. The crucial point in the 1 /N expansion is that the resummed a propagator has the non-perturbative 
behavior [25[ at large \p\, thus damping the degree of divergences with respect to usual perturbative graphs 
(such as those in Fig. 14) and thus the absence of momentum dependent divergences in Eq. (|C12|) . Nevertheless, 
since our construction relies on standard perturbation, in principle we would have to introduce corresponding new 
counterterms, with higher derivatives, into the Lagrangian at this perturbative order. But actually these p-dependcnt 
pieces contributions are canceled in the free energy graphs, as we will see next, so that it is not needed to deal with 
these extra non-renormalizable terms. Inserting the above two-loop counterterms into the (one-loop) free energy 
graph, for instance for the A counterterm in Eq. (|C9[) yields 



Kff _ , f dd p _ i 
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which gives, after taking the trace and going to Euclidean momenta 
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where we simply used for the integrand p 2 /(p 2 + ff) = 1 — rj 2 / (p 2 + tj 2 ). Note also that the integral J d n pp 2 is zero in 
dimensional regularization, so that only a mass dependent divergence remained in Eq. (|C14|) and we finally obtain 
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Similarly, for 



one obtains 



Note that the final mass and wave function counterterms are proportional to \ 2 rp , which simply reflects that the coupling A of the 
3-dimensional GN model has mass dimension —1 , i.e. the counterterms coefficients are actually dimensionless. 
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Next, the contributions in Eqs. (|C15|) , (|C17[) should be added in the MS scheme to the bare three loop free energy 
expression (IB4|) . This leaves a remaining divergent contribution: 
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(C18) 



which is finally rcnormalized by an additional vacuum (free energy) counterterm. This additional counterterm is 
expected, since the free energy is a composite operator, and is needed similarly in the (renormalizable) d = 1 + 1 GN 
model, as shown in a previous reference 13]. Note that, although the perturbative non-renormalizability manifests 
itself at order <5 2 by the presence of the higher derivative divergences in Eqs. ()C9|) . (|Cllj) . corresponding counterterms 
are not needed for the quantities (and perturbative order) we restrict ourselves to. Thus the arbitrariness in the final 
renormalized free energy is not more than the usual renormalization scale introduced from dimensional regularization. 
The complete renormalized three-loop contribution to the free energy is obtained, in the MS scheme, by adding all the 
finite contributions as given in Eqs. (|6.6[) . (|6.7p . (|C15|) . (|C17|) . respectively, leads to the final expression, Eq. 
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